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High-fidelity, efficient quantum nondemolition readout of quantum bits is integral to the goal of 
quantum computation. As superconducting circuits approach the requirements of scalable, universal 
fault tolerance, qubit readout must also meet the demand of simplicity to scale with growing system 
size. Here we propose a fast, high-fidelity, scalable measurement scheme based on the state-selective 
ring-up of a cavity followed by photodetection with the recently introduced Josephson photomulti¬ 
plier (JPM), a current-biased Josephson junction. This scheme maps qubit state information to the 
binary digital output of the JPM, circumventing the need for room-temperature heterodyne detec¬ 
tion and offering the possibility of a cryogenic interface to superconducting digital control circuitry. 
Numerics show that measurement contrast in excess of 95% is achievable in a measurement time of 
140 ns. We discuss perspectives to scale this scheme to enable readout of multiple qubit channels 
with a single JPM. 


I. Introduction 

Over the past decade, circuit quantum electrodynam¬ 
ics (cQED) has emerged as a powerful paradigm for scal¬ 
able quantum information processing in the solid state 
[1-4]. Here a superconducting qubit plays the role of 
an artificial atom, and a thin-film coplanar waveguide or 
bulk cavity resonator is used to realize a bosonic mode 
with strong coupling to the atom. Interaction between 
the qubit and the cavity is described by the Jaynes- 
Gummings Hamiltonian [5]. Strong interaction between 
the qubit and the cavity has been used to realize high- 
fidelity multi-qubit gates [6-9]; moreover, the qubit has 
been used to prepare highly nonclassical states of the 
resonator [10, 11]. In the limit where the qubit is far 
detuned from the cavity resonance so that A = wc — ojq 
satisfies |A| 3> gq, where ujc is the cavity frequency, ujq 
is the qubit frequency, and gq is the qubit-cavity cou¬ 
pling strength, the following dispersive approximation to 
the Jaynes-Gummings Hamiltonian is realized [1] (with 

;i= 1): 

= (wc+XQO'z)aia-i(a;Q - xq)o-^; (1) 

here XQ = 9q/^ i® dispersive coupling strength of 
the resonator to the qubit, and is the Pauli -2 opera¬ 
tor. One sees from the first term that the effective cav¬ 
ity frequency acquires a shift that depends on the qubit 
state. It is therefore possible to perform a quantum non¬ 
demolition measurement of the qubit by monitoring the 
microwave transmission across the cavity at a frequency 
close to the cavity resonance, for example, by using stan¬ 
dard homodyne or heterodyne techniques [1, 12]. This 
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approach for reading out the qubit state through cavity 
transmission measurements has become standard prac¬ 
tice. 

Recently much effort has been devoted to the devel¬ 
opment of near quantum-limited superconducting ampli¬ 
fiers for single-shot detection of the qubit state. Specific 
milestones include observation of quantum jumps in a 
transmon qubit [13], heralded state preparation of sin¬ 
gle qubit states to eliminate initialization errors [14, 15], 
deterministic preparation of entangled states [16], stabi¬ 
lization of qubit Rabi oscillations using quantum feed¬ 
back [17], and quantum teleportation [18]. The technol¬ 
ogy allows high readout speed [19] and entanglement over 
large distances [20]. While this approach works well for a 
small number of readout channels, the required supercon¬ 
ducting amplifiers, cryogenic semiconducting postampli¬ 
fiers, and quadrature mixers entail significant experimen¬ 
tal overhead: the amplifiers often require biasing with a 
strong auxiliary microwave pump tone which must be 
isolated from the qubit circuit with bulky cryogenic iso¬ 
lators; moreover, there is no clear path to integrating 
heterodyne detection at low temperature to provide for 
a compact, scalable architecture. 

An alternative approach that has not yet been consid¬ 
ered is to measure the state of the qubit using a photon 
counter. In contrast to an amplifier, which performs a 
linear mapping of input modes a, to output modes 
6,6/ a photon counter responds to the total power of 
the input signal a^a in a nonlinear fashion: the presence 
or absence of photons projects the counter into one of 
two possible classical output states, irrespective of the 
phase of the input signal. In the optical frequency range, 
the prototypical photon counter is the avalanche photo¬ 
diode [21, 22]: here, absorption of a single photon creates 
an electron-hole pair; the reverse bias of the pn junction 
sweeps the charge away from the depletion region and im¬ 
pact ionization generates additional electron-hole pairs. 
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FIG. 1: (Color online) (a) Schematic diagram of the Joseph- 
son photomultiplier (JPM) circuit. The junction is biased 
with a dc current, and microwaves are coupled to the junc¬ 
tion via an on-chip capacitor. In the simplest implementation, 
switching of the junction creates a voltage pulse that is read 
out by a room temperature comparator circuit, (b) Junction 
potential energy landscape. The junction is initialized in the 
ground state \g). An incident photon induces a transition to 
the first excited state | e), which rapidly tunnels to the contin¬ 
uum with rate yj. (c) Counter-based measurement in cQED. 
“Bright” and “dark” cavity pointer states result in binary digi¬ 
tal output from the JPM: “click” or “no click”, (d) In-phase (I) 
and quadrature (Q) phase space portrait of the cavity state 
after the ring-up, highlighting pointers to the |0) state in red 
and to the |1) state in blue. 


leading to a large and easily measured classical current. 

We have recently introduced a superconducting device 
that performs as a microwave-frequency analog of the 
avalanche photodiode [23-25]. The detector is a Joseph- 
son junction that is biased with a current such that the 
energy separation between the ground \g) and first ex¬ 
cited |e) states in the metastable minima of the junction 
potential energy landscape is resonant with the energy of 
the incident microwaves (see Fig. la-b). Absorption of a 
single microwave photon promotes the junction from the 
|g) to the |e) state, which tunnels rapidly to the contin¬ 
uum, producing a large and easily measured voltage of 
order twice the superconducting gap voltage. We refer 
to the detector as the Josephson photomultiplier (JPM). 
The JPM provides an intrinsically broadband frequency 
response; as we will show here, single-shot measurement 
contrast around 95% - suitable for scalable surface codes 
[26] - is achievable; the detector requires no microwave 
biasing, facilitating wireup of complex multi-qubit cir¬ 
cuits comprising many measurement channels; finally, 
the detector produces a binary digital output that in¬ 
terfaces well to scalable cold control circuitry based on 
Single Flux Quantum (SFQ) digital logic [27]. 

This paper is organized as follows. In Section II, we 
describe the basic principles of the JPM and discuss de¬ 
tector operation. In Section III, we present a detailed 
theoretical model of the proposed measurement protocol, 
with a focus on measurement contrast and back action. 
In Section IV, we discuss how close this scheme comes to 
a quantum non-demolition (QND) measurement, and in 


Section V we consider interactions with the environment, 
taking into account the full Jaynes-Cummings Hamilto¬ 
nian between the cavity and the qubit. Section VI is de¬ 
voted to a discussion of issues related to scaling this mea¬ 
surement approach to a large number of readout chan¬ 
nels. In Section VII we present our conclusions. 


II. Microwave Photon Counter Based on a 
Josephson Junction 

A schematic diagram of the JPM is shown in Fig. 
la. The Josephson junction is biased in the supercurrent 
state with a current If, that is slightly below the junction 
critical current Iq. The potential energy landscape U{5) 
for the phase difference <5 across the junction takes on 
a tilted-washboard form [28], with local potential min¬ 
ima characterized by a barrier height AU and plasma 
frequency ujp (Fig. lb). The circuit design and bias pa¬ 
rameters are chosen so that there is a handful of discrete 
energy levels in each local minimum of the potential; the 
JPM initially occupies the ground state \g). Microwaves 
that are tuned to the junction resonance induce a transi¬ 
tion to the first excited state |e), which rapidly tunnels to 
the continuum. This tunnelling transition in turn leads 
to the appearance of a large voltage across the junction 
of order twice the superconducting gap. Absorption of a 
photon thus yields an unambiguous and easily measured 
“click”. 

The experimental protocol involves pulsing the bias 
point of the JPM for a finite interval of order 10s of ns 
so that the transition frequency between the \g) and je) 
states is close to the frequency of the incident photons: at 
this point, the junction is in the “active” state, and there 
is high probability that absorption of a photon will induce 
a transition to the continuum. In the absence of resonant 
photons, there is a small, nonzero probability that the 
JPM will transition due to quantum tunnelling from \g), 
a dark-count event. JPM intrinsic contrast peaks for a 
bias such that AU/hiUp ~ 2 for a measurement interval 
that is roughly equal to the Rabi period of the coherent 
drive [25, 29]; for very short times, the interaction with 
the drive field is too weak to induce a transition, while for 
longer measurement times dark counts due to quantum 
tunnelling from the ground state degrade performance. 
In prior work, we have demonstrated efficiencies of order 
90% for coherent drive corresponding to Rabi frequencies 
around 100 MHz for junctions with extremely modest 
coherence times of order a few ns [23]. 

In the context of qubit measurement, the utility of the 
JPM is its ability to map bright and dark cavity states to 
two distinct classical output states: “click” or “no click”. 
It hence presents a measurement paradigm different from 
that of a linear amplifier and should be discussed in dif¬ 
ferent terminology [30]. For example, the gain of a JPM 
at an infinitesimal input signal is negligible as such a sig¬ 
nal will not activate it into the voltage state, whereas 
above a certain threshold the nonlinear gain is extremely 
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high. A performance comparison can, however, be done 
on the level of the overall qubit measurement protocol. 

In a conventional cQED measurement, the state of the 
qubit is encoded in the quadrature amplitudes of a weak 
microwave signal that is transmitted across the readout 
cavity. It is possible to access these amplitudes by pream¬ 
plifying the signal using a low-noise linear amplifier fol¬ 
lowed by homodyne or heterodyne detection; assignment 
of the detected signal to the qubit |0) or |1) states is per¬ 
formed by subsequent post-processing and thresholding. 
In the following, we analyze an alternative protocol in 
which the state of the qubit is mapped to the photon 
occupation of the cavity. The JPM then provides a high- 
fidelity digital detector of cavity occupation (see Fig. Ic- 
d). The measurement provides no information about the 
phase of the transmitted microwaves, or indeed about 
the amplitude of the transmitted signal beyond the dig¬ 
ital “click” / “no click” output of the JPM. As we show 
below, measurement contrast achievable with the JPM is 
comparable to that achieved with quantum-limited linear 
amplifiers, while the JPM provides unique advantages in 
terms of scaling to a large number of measurement chan¬ 
nels. We note that related proposals for photon counters 
were put forth recently that include both irreversible pho¬ 
ton absorption [31-33] and non-destructive photon de¬ 
tection via nonlinearity of a transmission line coupled to 
transmons [34, 35[. 
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FIG. 2: (Color online) Schematic of the three-stage measure¬ 
ment protocol: the upper panel describes the relevant con¬ 
trol, while the lower panel represents the corresponding cav¬ 
ity state, (a) In the drive stage, the cavity is driven strongly 
and coherently at the cavity frequency dressed by the qubit 
|1) state, oJd = cui = aic — XQ + XJ> for a duration td = tt/xq 
(assuming a square pulse). This projects the qubit onto ei¬ 
ther |0) or |1) and conditionally populates the cavity with a 
large number of photons n ~ [oiP when the qubit is pro¬ 
jected onto the |1) state, (b) During the measurement phase, 
the JPM is tuned into resonance with the cavity and allowed 
to interact; a bright cavity switches the JPM to its voltage 
state while a dark cavity leaves the JPM in the supercurrent 
state. This conditionally squeezes the cavity state by a small 
amount (not shown here), (c) In the reset stage, the cavity 
is again driven coherently at a;d, conditionally displacing the 
cavity to a near-vacuum state. 


III. cQED Measurement with a Microwave 
Photon Counter 

The basic scheme for qubit measurement with the 
JPM is shown in Fig. 2. The qubit (resonating around 5 
GHz) is coupled to a readout cavity (resonating around 
6 GHz). As in the usual dispersive limit of the Jaynes- 
Gummings Hamiltonian (1), the cavity acquires a disper¬ 
sive shift XQ = 9q/^ that depends on the state of the 
qubit. For the purposes of realizing a fast measurement, 
it is desirable to engineer a dispersive shift xq/tt « 10 
MHz, as opposed to the smaller dispersive shifts of or¬ 
der 1 MHz realized in typical cQED experiments. The 
measurement proceeds in three stages: (1) Drive stage. 
Here, we map the qubit state to microwave photon occu¬ 
pation of the readout cavity. A microwave pulse applied 
to the dressed frequency corresponding to qubit state |1) 
creates a “bright” cavity if and only if the qubit is in the 
excited state. If the qubit is in the ground state, the cav¬ 
ity acquires a non-negligible occupation at the start of 
the pulse, but it coherently oscillates back to the “dark” 
vacuum state upon completion of the drive pulse. Dur¬ 
ing the drive stage the JPM idles at a frequency that 
is blue detuned from the cavity by around 1 GHz. (2) 
Measurement stage. Here, we map photon occupation of 
the cavity to the voltage state of the JPM (“click” or “no 
click”). The JPM is rapidly tuned into resonance with 
the cavity. A bright cavity induces a transition to the 
voltage state, while a dark cavity leaves the JPM in the 


supercurrent state. (3) Reset stage. It is advantageous 
to coherently depopulate the bright cavity in order to 
circumvent the need for the cavity to decay via spon¬ 
taneous emission. However, since the depletion of the 
cavity due to interaction with the JPM is a stochastic 
process, so that neither the number of photons removed 
nor the back action on the cavity is perfectly known or 
reproducible, it is not possible to return the cavity pre¬ 
cisely to the vacuum state. Nevertheless, an appropriate 
coherent pulse can return the cavity to a state that is 
close to the vacuum. The measurement pulse sequence is 
shown in the upper panel of Fig. 2. 

In the dispersive regime of the qubit-resonator sytem, 
the unitary evolution of the full system is described by 
the Hamiltonian 

H = HeS + A{t) (a -I- -I- gj (aa^ + a'l'dj) , 

( 2 ) 

where u’j(t) is the frequency of the JPM, A{t) is the clas¬ 
sical drive applied to the cavity, gj is the cavity-JPM 
coupling, and iJeff is defined in Eq. (1). The JPM oper¬ 
ators couple the ground and excited state of the JPM, 
which are separated by a frequency wj (f) but do not cou¬ 
ple to the measured state. The JPM self-Hamiltonian 
contains df = diag(l, —1, F). Here, the energy of the 
measured state E is irrelevant once the tunnelling rates 
(which are not contained in this Hamiltonian as they re¬ 
quire interaction with an environment) have been fixed 
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independently. The measured state plays no role in the 
unitary dynamics of the system as it only couples inco¬ 
herently to all other states, and the full dynamics of the 
JPM are described by a Lindblad-type master equation. 

In the following we analyze the three stages of the mea¬ 
surement in detail. 


A. Drive Stage 

The goal of this stage is to prepare a photonic state 
in the cavity that is dependent on the qubit state, such 
that the conditional cavity states can later be distin¬ 
guished by the JPM in the measurement stage. The 
JPM idles in this stage, biased far off-resonance from 
the cavity such that the effective interaction between the 
cavity and the JPM is dispersive, with a dispersive shift 
Xj = ffj/(wc - wj). 

The effective Hamiltonian for the cavity becomes 

He = ojca^a + A{t) {a + a'^) , (3) 

where Cjq = (wc i XQ + Xj)- We choose a classical drive 
A{t) = ag cos (u>dt) for 0 < t < id where ag is the drive 
strength, wj the drive frequency, and tg the pulse length 
(for simplicity here we assume a square pulse). By setting 
Wd = oje ~ XQ + Xj obtain an effective cavity-drive 
detuning 6oj = cue ~ that depends on the state of the 
qubit: 


J 2xq for qubit in state |0) 

1 0 for qubit in state |1). 



FIG. 3: (Color online) Cavity occupation as a function of the 
duration of the applied drive, tg, for the qubit in states |0) 
and |1). Here xq/^ = 10 MHz, such that the optimal drive 
time is tg = 100 ns. 


where the qubit state |i/;)q = a|0) -I-&|1) can be prepared 
independently by the qubit drive line. After the drive 
stage, the system is left in the state 

l'l'(^d)) = (a|ao)c C) |0)q J- 6|q!i)c 0 |1)q) 0 |0)j, (7) 

which can be verified by solving Eq. (2) analytically. In 
the case that ag = 0, the cavity has nonzero occupation 
only when the qubit is in the excited state. 

The drive stage can be thought of as the first step in 
a quantum measurement of the qubit state, as described 
in the pointer basis formalism of Zurek [36]. In this lan¬ 
guage, the cavity states \ag/i) form the pointer basis that 
is entangled with the qubit. Examining the reduced den¬ 
sity matrix of the qubit state 


For such a classical drive of duration tg, the solution 
to Eq. (3) is easily obtained. Depending on the state of 
the qubit, the cavity will be in the coherent state |ao/i)i 
with 


ao = - l) ; ai =-'^tg (5) 

4xq 2 

up to a global phase. We see that when the qubit is in 
state |0), the cavity occupation oscillates sinusoidally at 
a frequency set by the detuning of the drive pulse from 
the dressed cavity resonance. On the other hand, when 
the qubit is in state |1), the cavity occupation |aip grows 
monotonically in time. This is shown in Fig. 3, where 
we plot cavity occupation versus coherent drive time for 
the qubit in states |0) and |1). In order to maximize con¬ 
trast between the dark and bright cavity states to which 
the qubit states |0) and |1) are mapped, it is optimal to 
choose tg = tt/xq such that ao = 0 at the end of the 
drive stage. The length of the drive stage is therefore set 
by the requirement that ag{tg) = 0 and not by the input 
cavity coupling, which is the inverse of the decay time of 
the cavity through its input port. 

We assume the system starts in the state 


( 6 ) 


. _ f |a|2 a*bD \ 
y ab*D |&|2 ) ’ 


( 8 ) 


we see that qubit coherence has been suppressed by a 
factor D = exp (—|ai — aoP)j which quantifies the de¬ 
phasing of the qubit induced by the interaction with the 
pointer basis (cavity). The dephasing would be complete 
if the pointer states were orthogonal. Moreover, mapping 
of the qubit |1)/|0) states to bright/dark cavity states can 
be viewed as a coherent amplification step, as the infor¬ 
mation about the qubit state is now contained in a large 
number of photons. A more detailed discussion of the 
consequences of this overlap on the detection contrast 
and back action will be presented later. 

As a result of the strong dephasing of the qubit state 
during the drive stage (quantified by the factor D), our 
multi-stage protocol explicitly exposes the role of the pre¬ 
measurement stage in quantum non-demolition (QND) 
readout. In particular, our protocol highlights the fact 
that in QND readout of the qubit state, measurement 
of the cavity pointer states is not the major source of 
qubit state dephasing. The qubit states are dephased 
during the pre-measurement, when qubit states and cav¬ 
ity pointer states are entangled, which in our case corre¬ 
sponds to the drive stage. The main role of the subse¬ 
quent pointer state measurement (the measurement stage 


14/(0)) = |0)c®|V')q®|0)j, 
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in our protocol) is to break unitarity and “freeze” the 
qubit in a dephased state. This distinction between pre¬ 
measurement and measurement is less obvious in qubit 
readout using a continuous cavity signal with linear am¬ 
plification and heterodyne detection. The clear distinc¬ 
tion between pre-measurement and measurement in our 
protocol allows for independent control of each stage, 
which can be used to achieve higher readout fidelity (as 
has been done here), and to study, both in theory and 
experiment, QND measurement and the pointer basis for¬ 
malism with an explicit physical system in mind. A sim¬ 
ilar distinction between pre-measurement and measure¬ 
ment exists in a readout scheme in atomic cavity QED, al¬ 
beit in a rather different parameter regime. This scheme 
employs dispersive coupling between the cavity and a 
travelling atom (pre-measurement) followed by atomic 
state detection via ionization (measurement) to read out 
the cavity state [37, 38]. 


B. Measurement Stage 

After the drive stage, the qubit state information has 
been transferred to the cavity occupation. In the mea¬ 
surement stage, a measurement of the cavity by the JPM 
will reveal the state of the qubit. During this stage, 
the JPM is brought into resonance with the dressed fre¬ 
quency of the cavity corresponding to the qubit |1) state, 
(jjj = (jjc — XQ, in order to maximize detection in the 
case that the qubit is excited. In practice, precise tuning 
of the JPM bias point is not required due to the broad 
detection bandwidth of the JPM [25]. 

The Hamiltonian during this stage is that of Eq. (2) 
with A(t) = 0. In the following, we assume a cavity- 
JPM coupling = 50 MHz. In addition, the system 

evolves incoherently as a result of tunnelling (both bright 
and dark) and relaxation of the JPM. We consider tun¬ 
nelling from both the JPM excited and ground states 
to the measured state, and relaxation from the excited 
state to the ground state, with corresponding rates yj, 
7d, and 7 r, respectively. Here we take yj = 200 MHz, yo 
= 1 MHz, and yp = 200 MHz; this relaxation rate cor¬ 
responds to a junction with capacitance 100 pF directly 
connected to an environmental impedance of 50 D. The 
total evolution of the system can therefore be described 
by a Lindblad-type master equation with Lindblad op¬ 
erators corresponding to each incoherent process of the 
JPM, as outlined in more detail in our previous work 
[24, 25]. 

As the cavity-JPM coupling and bright count rate can 
be controlled independently of one another, they can be 
adjusted into an optimal regime for good measurement. 
As explained in more detail in our previous work [23, 
25], the optimal regime for good measurement is when 
5 j ~ 7J) in this regime the bright count rate is large 
enough for a bright count to occur within the occupation 
time of the JPM (per Rabi cycle), while not so large 
as to result in a Zeno effect suppression of the cavity- 




FIG. 4: (Color online) (a) Detection probability P (|ap,tni) 
as a function of measurement time tm (horizontal axis) and 
cavity occupation |ap (color) for system parameters gj/2Tv = 
50 MHz, yj = 200 MHz, yo = 1 MHz and yn = 200 MHz. 
(b) Detection probability P (|ap,tni) as a function of cavity 
occupation |ap, both numerical simulations and analytic fit, 
for tm = 25 and 50 ns. 

JPM interaction. The coupling and bright count rate 
chosen for the numerical simulations presented here are 
well within the optimal regime for good measurement. 

Starting from the output state of the drive stage, 
we numerically solve the master equation for the mea¬ 
surement stage to obtain the detection probability 
P (jap, tin) as a function of cavity occupation and mea¬ 
surement time tin- In the case that the qubit starts in 
the |0)/|1) state, the detection probability reduces to 
P (|ao/iP,tm)- We define the qubit measurement con¬ 
trast as the difference in detection probability between 
these two cases: 

C = P(|ai|2,t„,)-P(|ao|^^m). (9) 

Clearly, the measurement contrast is optimized when 
P (jaop, tm) = 0, which requires that ao = 0 and yo = 0. 
The measurement contrast has a maximal value of one if 
jaij —)■ oo, indicating a perfect measurement. 

Figure 4(a) shows the detection probability as a func¬ 
tion of both jap and tm for the system parameters pre¬ 
viously discussed. For all coherent states with average 
photon number jap > 0 we see similar behavior as a 
function of time, with saturation of the detection prob¬ 
ability around 40 ns, irrespective of cavity occupation. 
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The fact that the detection probability saturates at a 
value less than unity is explained by the two competing 
mechanisms for excitation loss in the JPM: measurement 
tunnelling and inelastic relaxation, only the former of 
which results in a bright count. The black curve in Fig. 
4(a), a = 0, is the detection probability for the cavity 
when the qubit is in the ground state, P (|aoP = 
and the fact that it is nonzero is due only to dark counts, 
which occur with a probability PD(0,t) = 1 — 

In a simplified picture where energy in the detector 
automatically leads to a click, we would have 

P(|an = PedaH + (1 - PB(|an)PD, (10) 

where Pedap) is the bright count probability and Pd = 
Pd (0, t) is the dark count probability, which we take to be 
independent of a. However, the detection mechanism of a 
JPM involves the coherent absorption of energy prior to a 
tunnelling transition to the classically observable voltage 
state. In the interval following absorption of a photon 
but prior to tunnelling, dark counts cannot occur as the 
JPM is in its excited state. This breaks the dark/bright 
symmetry of Eq. (10); as a result, this equation may no 
longer be valid. However, as it shows how dark counts are 
detrimental to measurement contrast, Eq. (10) is a good 
reference point to compare against, and will still be valid 
in some situations, such as when the JPM coupling rate 
is smaller than the bright tunnelling rate, i.e., gj < jj. 

Figure 4(b) shows the detection probability as a func¬ 
tion of |ap for <ni = 25 and 50 ns, along with an analytic 
fit to the data by Eq. (10), with Pb given by the curve 

PB(|Q!|^,tm -t oo) = 1 - exp (-\af —^ . (11) 

V 7J + 7R / 

See Appendix A for a derivation. The analytic curves for 
both = 25 and 50 ns are so similar on this scale that 
only tni = 50 ns is plotted. As can be seen, the analytic 
fit is valid when is sufficiently large. For small tm, Eq. 
(10) remains close to correct, but the approximation for 
Pb in Eq. (11) breaks down. 

We have calculated detection probability P(|ap,tin) 
and measurement contrast versus measurement time tm 
for |aoP = 0, and |aip = 10; the results are shown in 
Fig. 5(a). The measurement contrast peaks at « 95% 
around 40 ns, indicating that a good choice for is 40 ns. 
At longer times the measurement contrast will eventually 
begin to decrease, as P (laoPjfm) continues to increase 
while P (|aip,tni) asymptotes to near unity. In general, 
we observe that increasing ai increases the contrast, ul¬ 
timately limited by the breakdown of the dispersive ap¬ 
proximation to the Jaynes-Cummings Hamiltonian. 

The contrast shown in Fig. 5(a) is for one set of sys¬ 
tem parameters, and in principle it is possible to obtain 
higher values of measurement contrast by optimizing over 
parameter space. Figure 5(b) shows the optimal mea¬ 
surement contrast as a function of bright count rate, yj, 
for various bright states |aip. The ratio of the bright 
and dark count rates is set by fabrication parameters of 




FIG. 5: (Color online) (a) Excited qubit detection prob¬ 
ability P(|aip,tm), ground qubit detection probability 
P (|ao|^, tm), and measurement contrast versus measurement 
time. Here |ai|^ = 10 and |ao|^ = 0, with system param¬ 
eters as before, (b) Measurement contrast as a function of 
bright count rate yj, for various |ai|^ and = 50 ns. The 
relaxation rate remains fixed at yn = 200 MHz, while the 
dark count rate changes such that the ratio yj/yo = 200 is 
unchanged. 


the JPM, and therefore remains fixed at yj/yD = 200. 
However, the inelastic relaxation rate remains fixed as yj 
changes, such that y^ = 200 MHz regardless of the value 
of yj. As can be seen in Fig. 5(b), within experimentally 
reachable parameter regimes contrast greater than 95% 
is possible. 

Measurement contrast is ultimately limited by the pos¬ 
sibility of misidentifying the qubit state. Misidentifica- 
tion of the excited state as the ground state is due to the 
nonzero vacuum component of the coherent state jai) as 
well as to internal photon loss. This occurs with a proba¬ 
bility 1 —P (|aiP, tm), which is bounded below by 
(occuring when yn = 0). Misidentification of the ground 
state as the excited state is the result of a dark count (as¬ 
suming ao = 0), and the probability of misidentification 
in this case is exactly PD(0,tm) discussed earlier. The 
problem of misidentification, and the fact that measure¬ 
ment contrast is less than unity even for yR,yD = 0, is 
related to the basis of our measurement protocol and will 
be discussed in more detail in Section IV. 

After the measurement stage, if the JPM absorbs a 
photon and switches out of the supercurrent state, clas- 
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sical emission due to this switching process could induce 
relaxation in the qubit or produce a spurious population 
in the readout cavity that would spoil the reset pulse 
[39]. The resulting population is proportional to the en¬ 
ergy spectral density of the classical current drive at the 
qubit or cavity frequency. A straightforward approach 
to address this would be to install a microwave isolator 
between the cavity and JPM, as in conventional cQED 
experiments, where one inserts one or more cryogenic iso¬ 
lators between the measurement cavity and the super¬ 
conducting preamplifier. The breaking of time-reversal 
symmetry by the isolator allows signals to travel from 
the cavity to the readout device with minimal loss, while 
back action noise is heavily attenuated. However, we 
anticipate that this classical back action can be greatly 
suppressed by an appropriate choice of JPM parameters 
to suppress harmonics of the switching transients at the 
qubit and cavity frequencies, or by shunting the JPM 
by an appropriate admittance to prevent a full switch of 
the JPM phase to the running state. Alternatively, it 
might also be possible to eliminate a cryogenic isolator 
by incorporating on the JPM chip tunable impedance¬ 
matching circuitry, as this would allow for the realiza¬ 
tion of a strong impedance mismatch between the cavity 
output and JPM immediately after the end of the mea¬ 
surement stage. 


C. Reset Stage 

The final stage is to remove the energy from the cavity, 
ideally leaving the cavity-qubit system in the conditional 
states |0)c|0)q or |0)c|1)q to allow for additional opera¬ 
tions on the qubit. This can be achieved through cavity 
decay by simply waiting long enough; however, because 
the total cavity decay time may be comparable to the 
qubit Ti, it is preferable to actively reset the cavity. 

After the measurement stage^ the cavity is either in the 
vacuum state or the state Here, Hjp are 

the back action operators [24] on the cavity due to JPM 
tunnelling and inelastic relaxation, respectively, and N 
is the number of photons removed from the cavity by the 
JPM. These back action operators interpolate between 
the standard lowering operator B = a and the subtrac¬ 
tion operator 13 = arT^!'^ [24]. We neglect for the mo¬ 
ment the classical back action on the cavity due to the 
transient current that develops when the JPM switches 
to the voltage state. As a starting point for reset, we 
will assume that 13^13^~^\ai) ~ |a]y[) even with large 
7 r,, with 


|q;m|^ = Tr 




( 12 ) 


the average photon number of the cavity state after mea¬ 
surement. 

At the end of the reset stage, we desire the cavity to 
be in the vacuum state independent of the qubit state, 
and thus we must invert the drive stage. Consider a 


Hamiltonian of the form of Eq. (3), with a more general 
drive A{t) = ai cos (wdt + (/>)0(td — t), where td = tt/xq 
as before. The unitary operation applied to the cavity is 
then 


C, =IcO|0)(0|q + .D(/3)o|1)(1|q. (13) 

Here D{/3) is the displacement operator on the cavity, 
with 



Thus, by choosing ai such that (aitd)/2 = |q;m| and 
setting (j) = (2n + l)7r,n G Z, we have (3 = —au- Under 
these conditions, the operation 14 will leave the cavity 
in the vacuum state independent of the qubit state, and 
will do so with an operation time A = td, significantly 
shorter than the total decay time of the cavity. 

However, after detection by a JPM the state of the 
cavity is not a coherent state; thus there does not exist a 
displacement operator D{/3) that will map it identically 
to the vacuum state. One can quantify the resulting de¬ 
viation from vacuum by calculating 

Eiai,N) = 1 - ±\{0\D{-aM)BdB^-^\ai)\^ (14) 


where A is the normalization of the state after measure¬ 
ment. This error will depend on the form of the back 
action. Assuming all back actions can be expressed in 
terms of subtraction operators as in [24], we find 


E{ai,N) = 1- —\{-aM\B^\ai)\‘^ 


= 1 - 


P, 


N 


\^M I +1^11 


E 


n ^n+N 


^n\{n + N)\ 


(15) 


Here, the normalization A^ is the probability of N pho¬ 
tons being subtracted [40], Pn = 1 — where 

r(A^, jaip) is the upper incomplete Gamma function. 
The error of this reset pulse is shown in Fig. 6 for dif¬ 
ferent values of N and as a function of jaip. As can be 
seen, the maximal error increases with increasing N, but 
for all N the error tends to zero as jaij —)■ oo. In reality, 
as the value of N is not fixed, a better estimate for the 
average error can be obtained by averaging over the er¬ 
ror traces shown in Fig. 6. Note that if the back action 
operator is closer to the standard photon lowering oper¬ 
ator a, this figure of merit will improve as the coherent 
states are eigenstates of this operator, and can be moved 
to vacuum exactly. 

The possibility exists that more complicated pulse se¬ 
quences during the reset phase will be able to map the 
cavity state identically to the vacuum; however, consider¬ 
ation of such sequences is beyond the scope of this work. 
In any case, the error of the reset pulse shown here does 
not affect the success of qubit readout. 
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function of the measurement time tm) using the following 
Choi matrix fidelity 



FIG. 6: (Color online) Numerically evaluated overlap error for 
the reset pulse described in the text. The number of photons 
N removed from the cavity is shown in the legend. 


IV. QNDness of the Measurement 


Ideally, we would like our measurement to project the 
qubit into its bare basis, {|0/1 )q}, hence implementing 
a quantum nondemolition (QND) measurement. A hall¬ 
mark of QND measurement is that a repeated measure¬ 
ment leads to the same result with certainty. Our mea¬ 
surement scheme starts from the dispersive Hamiltonian 
of Eq. (1) in the cavity ring-up phase, which is QND in 
the sense that the qubit and the pointer coupling com¬ 
mute. However, the process of destructive photon ab¬ 
sorption necessarily results in a deviation from QNDness, 
which we analyze in detail below. Even in the case of an 
ideal measurement (yu, 7 d = 0 and tm —> cx)), the poten¬ 
tial misidentification of the two states leads to a QND 
error. Starting from the state in Eq. (7), the measure¬ 
ment projects the qubit conditionally onto the states 


I,/, 1 - «e -f 6e 

“ y'|a|2e-l“ol= + |6|2e-l“d^‘^ ‘ 

^ _ aVl - e~l“ol^|0)Q -I- bVl - e~l°^il^|l)Q 
“ v/|a|2(l-e-l“olQ + |&|2(i_e-DilQ ' 


(16) 

(17) 


Even for uq = 0 these states are non-orthogonal (and not 
equal to the ideal QND post-measurement states), and 
their overlap is related to the overlap of the cavity pointer 
states luo/i)- This overlap is what allows misidentifica¬ 
tion to occur, ultimately limiting the contrast and QND¬ 
ness. 

If we do not condition on the measurement outcome, 
the effect of a perfect QND measurement is the quantum 
process defined by the map 

a|0)+6|l)^|an0)(0| + |6ni)(l|, (18) 


which completely destroys all coherence in the qubit 
state, while maintaining relative populations. We can 
describe this quantum process by its Choi matrix Cper 
(see Appendix B), and can compare this to the Choi ma¬ 
trix Ct^ describing our measurement protocol (which is a 



(19) 


As this fidelity compares the unconditional measurement 
protocol, it does not contain information about the suc¬ 
cess of the measurement (which we believe is well de¬ 
scribed by the contrast), but instead quantifies how close 
the possible qubit output states are to the desired ones. 
As a result, by choosing ideal QND measurement as our 
reference process, we can directly quantify the QNDness 
of our measurement protocol. 

As we are examining the unconditional measurement 
process, any measurement time dependence in will be 
due to changes in the back action of JPM measurement 
on the cavity that change the post-measurement cavity 
state, and therefore modify the coherence of the post¬ 
measurement qubit state. However, for ai ao suffi¬ 
cient decoherence of the qubit state has occurred during 
the drive stage that the resultant measurement time de¬ 
pendence of the fidelity is several orders of magnitude 
smaller than the average value, and in practice we can 
set = J-ao- Using Eqs. (16) and (17) we can calcu¬ 
late J-ao analytically in the ideal case when 7 r,7d = 0 
(see Appendix B): 

.Too = ^ (l + \/l-^(ao,ai)2) (20) 

K{ao,ai) = e-Kl“ol"+l“ilQ 

+ .y(l_ e-DoF) (i_e-DiF). 

For |aoP = 0 and |aip = 4, we already have J'oo > 
99%. For jofip = 10 as used elsewhere, Too > 99.99%. 
Ultimately, this value of the fidelity should be considered 
the fundamental limit of our protocol as it corresponds 
to the ideal case, ignoring both JPM relaxation and dark 
counts, as well as other environmental interactions. 

When JPM relaxation is non-negligible (yu yf 0), even 
for —>■ oo the measurement conditionally projects the 
qubit onto mixed states rather than the pure states of 
Eqs. (16) and (17), as even for ao = 0, Itpo) is mixed 
incoherently with a |1)(1|q component. Similarly, dark 
counts cause mixing of the state |^i) as they incoherently 
add a |0)(0| component to \ijji). Therefore, to describe 
back action, we can use POVM (positive-operator val¬ 
ued measure) elements for the qubit state to describe the 
map onto the post-measurement state. While full deter¬ 
mination of these POVM elements is beyond the scope 
of this work, the unconditional quantum process is 
also directly affected by changes in the POVM elements, 
such that it is quantitatively different when yu, 7 d yf 0. 
However, the average value of the fidelity is nearly the 
same, and as changes to the fidelity with measurement 
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time are several orders of magnitude smaller than the av¬ 
erage value, the fidelity is not a good measure to compare 
the ideal case with that for 7 r,7d 0. 

Therefore, to qualitatively study the deviations from 
QNDness introduced by JPM relaxation and dark counts, 
we examine the probability that repeated measurements 
(within qubit Ti) will give the same measurement result. 
Consider single measurement probabilities Pa, where a € 
{0,1} is the measurement outcome, and joint measure¬ 
ment probabilities Pab where a, b are the outcomes of the 
second and first measurements, respectively. For an ideal 
QND measurement as defined above, we have 

-Poo = pQ 

Poi = P ’10 = 0 

P ’11 = Pi- (21) 

When JPM relaxation and dark counts are taken into 
account, none of these relationships hold. This is gener¬ 
ally a result of the fact that our protocol can misidentify 
the qubit state (due to dark counts, energy relaxation, 
or less than unit contrast), so that the second event does 
not occur with unit probability. In particular, due to 
dark counts Pq ^ Pooj similarly, due to the probabil¬ 
ity of not detecting a photon for a given cavity state or 
mistakenly measuring the vacuum component of the |ai) 
state, we have Pn ^ Pi. Pn can also be further modified 
by imperfect reset of the cavity. The symmetry between 
Pqi and Piq is not broken by misidentification; however, 
they are both nonzero. On top of these misidentifications, 
QNDness can also be limited by corrections beyond the 
dispersive Hamiltonian as discussed in the next section. 

V. Environmental Interactions and Corrections 
Beyond the Dispersive Hamiltonian 

So far the discussion has focused on a closed qubit- 
cavity subsystem. When we consider interactions with 
the environment, it is apparent that the dominant effects 
are qubit and cavity relaxation. The timescale of these 
effects depends heavily on the frequency of the JPM, as 
it is most strongly coupled to environmental modes. 

During the drive and reset stages of the measurement 
protocol, the JPM idles at a frequency that is far blue 
detuned from the cavity or qubit resonances. As a re¬ 
sult, the leading order decay channel for both the cavity 
and the qubit is through the cavity’s input port. If we 
take a cavity decay rate k « 100 kHz, we find a qubit 
lifetime limitation of Tf « 2 ms for vacuum in the cav¬ 
ity [41] (see Appendix C for further details), and this Tf 
will in fact increase with higher cavity occupation [42], 
though this is a higher order effect not considered in our 
evaluation. This Purcell limited qubit Ti can only be 
calculated by first considering the full Jaynes-Cummings 
Hamiltonian when deriving the master equation for the 
coupled system, and so it is inherently not contained in 
the dispersive picture. Essentially, induced qubit Ti can 


be understood by observing that the eigenstates of the 
full JC Hamiltonian are always dressed (albeit weakly at 
strong detuning), and thus decay of the dressing cloud 
can lead to decay of the eigenstate. 

In addition, while spontaneous emission of the cavity 
through its input port is inconsequential, emission toward 
the JPM during the drive stage will degrade the prepara¬ 
tion of the cavity pointer states, with the dominant effect 
being a nonzero occupation japP of th® |0)Q-state pointer 
upon completion of the drive stage. However, this effect 
is very small due to the large cavity-JPM detuning dur¬ 
ing the drive stage, and so it only minimally affects the 
contrast. 

During the measurement stage, the JPM is brought on 
resonance with the cavity, and cavity decay through the 
JPM is desirable, since it amounts to bright tunnelling 
or JPM relaxation. However, as the cavity-JPM states 
hybridize, qubit decay through the JPM is also possible, 
as a result of beyond-dispersive effects identical to those 
for qubit Purcell decay discussed previously. Through a 
procedure similar to that of [41], we obtain a JPM-limited 
qubit lifetime of « 2 /xs during the measurement 
stage for vacuum in the cavity, considerably shorter than 
Ti « 2 ms (see Appendix C for further details). For 
an occupied cavity the situation is more complex, due 
to additional excitations as well as stimulated emission 
channels, but we find that to lowest order in gq/A the 
qubit lifetime increases as cavity occupation increases, 
and that for jap = 10 one would expect a qubit lifetime 
of « 40 ^s. 

However, numerical simulations (see Fig. 7(a)) in¬ 
dicate that there is no appreciable qubit decay proba¬ 
bility during the overall measurement process. We at¬ 
tribute this to the fact that the global state of the sys¬ 
tem is frozen once the JPM is in the measured state, 
and since 7 j ^ 1/T^^, this occurs long before any ap¬ 
preciable qubit decay. In fact, we expect only a 0.05% 
change in the qubit state due to JPM-mediated decay 
(see Appendix C for further details), which is completely 
washed out by other effects in Fig. 7(a). In other words, 
a seemingly short induced lifetime during the measure¬ 
ment stage is inconsequential if the associated relaxation 
channel is only open for a short time. This implies that 
a working point with a very fast bright tunnelling rate is 
optimal. 

The qubit also experiences dephasing due to low fre¬ 
quency noise at the JPM with a characteristic timescale 
However, as the ideal cQED measurement protocol 
should maximally dephase the qubit state, this low fre¬ 
quency noise does not affect the fidelity or measurement 
contrast of our protocol. 

To quantify measurement degradation due to beyond- 
dispersive effects, we compare the process fidelity of Eq. 
(19) for dispersive qubit-cavity coupling with that for the 
full Jaynes-Cummings Hamiltonian. Figure 7(b) shows 
the fidelity as a function of measurement time, for sim¬ 
ilar parameters as used throughout and joip = 9. As 
expected, the dispersive fidelity changes only minimally 
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Measurement Time (ns) Measurement Time (ns) 


FIG. 7: (Color online) (a) Qubit expectation for the qubit 
initially in the excited state versus measurement time (Jaynes- 
Cummings Hamiltonian) and (b) process fidelity of qubit 
readout for both the dispersive Hamiltonian and the Jaynes- 
Cummings Hamiltonian. Coupling strength and tunnelling 
rates are as used throughout. Here |aop = 0 and |ai|^ = 9. 

as a function of measurement time, while the Jaynes- 
Cummings fidelity both oscillates and grows with mea¬ 
surement time. Crucially, the fidelity for the full Jaynes- 
Cummings Hamiltonian is still approximately 98%, i.e., 
not significantly less than for the dispersive Hamiltonian. 
It is the focus of future study to improve this number. 

Finally, we have examined both the cavity occupation 
during the drive stage and the measurement contrast for 
the full Jaynes-Cummings Hamiltonian. For the drive 
stage the major effect is a small shift in the time td at 
which the cavity occupation is minimized for the qubit 
in the ground state, and an increase in the minimum 
occupation |aoP (Fig. 8(a)). This results in a reduction 
of the contrast (as can be seen Fig. 8(b)); however, this 
reduction in contrast is not significant enough to seriously 
degrade the success of our measurement protocol. 


VI. Scalability of Counting Measurement 

A useful multiqubit processor comprising hundreds if 
not thousands of qubits will require a large number of 
measurement channels with their associated wiring, fil¬ 
tering, and isolation. It is therefore important not only 
to examine the ultimate performance of a single measure¬ 


FIG. 8: (Color online) (a) Cavity occupation during the drive 
stage for the full Jaynes-Cummings Hamiltonian, (b) Mea¬ 
surement contrast for the dispersive Hamiltonian and the 
Jaynes-Cummings Hamiltonian. As elsewhere, xq/tt = 10 
MHz, srj/27r = 50 MHz, yj = 200 MHz, 7 d = 1 MHz and 
7 r = 200 MHz. In both plots, the drive strength is chosen 
such that |qiP = 9 for td = 100 ns for the dispersive Hamil¬ 
tonian. 


ment channel, but also to consider prospects for scaling 
to many measurement channels. From the standpoint 
of scalability, JPM-based counter measurement possesses 
several unique advantages compared to conventional het¬ 
erodyne measurement based on low-noise superconduct¬ 
ing amplifiers. 

The JPM requires only relatively low-bandwidth dc 
wiring for biasing, thus eliminating the need for bulkier 
cryogenic coaxial lines and microwave attenuators. More¬ 
over, operation of the JPM requires no microwave pump 
tone, eliminating a major source of cost, complexity, and 
deleterious crosstalk in conventional cQED circuits. In 
addition, because the output signal of the JPM is of the 
order of twice the superconducting gap, no cryogenic am¬ 
plifiers are needed and the JPM signal can be detected 
with straightforward room-temperature electronics. Al¬ 
ternatively, the binary digital output of the JPM pro¬ 
vides a natural interface to the SFQ-logic family. Here, 
classical bits are stored in the form of quantized volt¬ 
age pulses whose time integral equals the superconduct¬ 
ing flux quantum $0 = /i/2e. Optimized SFQ circuits 
can be clocked at 100s of GHz, and they offer orders 
of magnitude lower dissipation than conventional CMOS 
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logic. The integration of a classical SFQ control circuit 
in the multiqubit cryostat would yield significant reduc¬ 
tions in power consumption, latency, and overall system 
footprint. 

The large intrinsic bandwidth of the JPM (approaching 
1 GHz) [25] also allows for the possibility of time-domain 
multiplexing. Multiple qubits, each with a separate read¬ 
out cavity at slightly different frequencies, could be inter¬ 
rogated with a single JPM by selectively addressing each 
cavity with drive pulses at different frequencies. While 
the “click”/“no click” output of the JPM does not enable 
frequency-domain multiplexing of the cavity readout, it 
is possible to multiplex instead by staggering the read¬ 
out of the cavities in time, with an offset between cavity 
measurements of order 10s of ns. 


VII. Conclusion 

In conclusion, we have outlined a new readout scheme 
for superconducting quantum bits using selective cavity 
ring-up and photodetection. We show that even with¬ 
out detailed optimization, our measurement protocol is 
compatible with the requirements of fault tolerance, with 
achievable measurement contrast greater than 95% in 
measurement times of order 100 ns. Counter-based qubit 
measurement possesses distinct advantages in terms of 
scalability, with simple wireup and dc biasing require¬ 
ments and the prospect of multiplexing in the time do¬ 
main. Finally, as the counter maps quantum information 
to a binary digital output without the need for room- 
temperature heterodyne detection and post-processing, 
our scheme provides a natural interface between a su¬ 
perconducting quantum processor and cryogenic classical 
control circuitry based on the SFQ digital logic family. 
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Appendix A: Derivation of the Analytic Expression 
for the Bright Count Rate 

Following [25], we begin by assuming we have a dark 
count-free JPM coupled to a cavity in an A-photon Fock 
state. A single photon in the cavity would cause a bright 
count with probability Pi = 7j/(7r-I-7j), where ■jj is the 
bright tunnelling rate and 7 r is the inelastic relaxation 
rate of the JPM, as defined before. However, if instead 
the JPM relaxes back to the ground state, the second 


photon in the cavity will cause a bright count with prob¬ 
ability P 2 = [ 7 r/(7r+7j)]Pi, where the first factor is the 
probability that the first photon is lost due to inelastic 
relaxation. Therefore, for the nth photon, we have 


Pn 


/ \ "“1 
/ 7R \ 7J 

V7R + 7J/ 7r + 7j' 


Summing up all probabilities for n = 1,..., A, we obtain 

Pff = 1- f —\ = 1 - exp (-N In (^1 -k —^ ^ . 

V7R + 7J/ V V 7r// 


For a coherent state, we improve the estimate by averag¬ 
ing over A for a state with given jap: 


^(i«n = E 


|2Af 


A! 


-GIV 


N 


= 1 — exp — I a 


7J 


(7J + 7r) / 


(Al) 


This analytic expression is valid for 7 r > ■jj, assuming 
that all rates are independent of the number of photons 
in the cavity. This assumption implies that the photon 
excites the JPM faster than both rates 7 r and yj, and is 
valid for long measurement times. 


Appendix B: Analytic Derivation of Fidelity in the 
Ideal Case 


To calculate Poot we begin by defining the uncondi¬ 
tional map on the qubit state (with 7d,7r = 0) using 
Eqs. (16) and (17) by 

foo mm = Poimm+Piimm, (bi) 

for an arbitrary initial qubit state |^/>) = ajO) -I- 5|1). The 
“click”/“no click” probabilities for such an input state are 
given by 

Po = -k (B2) 

Pi = |a|2(l - e-l“»l') + |&|2(1 - e-l“il'). (B3) 

In light of this, Eq. (Bl) becomes 

£^mm = \am{o\ + \bm{M ( bj ) 

-k ab*K{ao, ai)|0)(l| -k a*bK{ao, q;i)| 1)(0|, 
where A(ao,ai) is as defined in Eq. (20): 

K{ao,ai) = e-HKP+UiP) + y'(i _ e-|a„k) (1 _ e-UiP). 

Now that the map is fully determined, we can calculate 
the Choi matrix elements 



ijkl 




(B5) 
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and find that the Choi matrix is given by 

1 0 0 K{aQ,ai) 

0 0 0 0 

0 0 0 0 

K{aoTai) 0 0 1 

By a similar procedure, the Choi matrix for perfect QND 
measurement is given by 




/ 1 0 0 0 \ 

I 0 0 0 0 I 

0 0 0 0 

\ 0 0 0 1 / 


(B7) 


With both Choi matrices defined, using Eq. (19) we can 
calculate the fidelity to be 

^oo = ^ (l + v/l-^(ao,ai)") , (B 8 ) 

as in Eq. (20). 


Appendix C: Cavity and JPM-Limited Qubit 
Lifetimes 

We first summarize the results of [41] for a qubit cou¬ 
pled dispersively to a cavity with decay rate k; next we 
extend this result to include the JPM. For a qubit coupled 
dispersively to a cavity, the dressed qubit-cavity eigen¬ 
states to second order are 

|l,n-l)« (l-^j |l,n-l)-^|0,n), (Cl) 

(l-^) |0,n) + ^|l,n-l), (C2) 

where | 0 / 1 , n) are the uncoupled eigenstates of the cavity- 
qubit system and A = ojq — ojq. We are interested 
in Purcell-limited qubit relaxation, i.e. transitions from 
|l,n) to | 0 ,n) mediated by the cavity’s coupling to the 
external environment, which we assume takes the stan¬ 
dard form with the cavity coupling operator given by 
X = a + . From [41], the decay rate for this process is 

given by 

rWo;i = ^(A__)|^X^P; (C3) 

here k(w) is the coupling constant that depends on 
the spectral density of the cavity’s environment, which 
should be evaluated at A^^ = + XQ- To lowest 

order in (/q/A, 

(C4) 

Assuming an Ohmic spectral density and using as a ref¬ 
erence value the coupling constant at the uncoupled cav¬ 
ity frequency k{ujc), we have k{uj) « k(wc)w/wci where 


LOc is the bare cavity frequency, not the ultraviolet cut¬ 
off frequency of the Ohmic spectral density [29]. Set¬ 
ting k(ojc) = 100 kHz (the input coupling to the cav¬ 
ity), = 1 GHz, and a corresponding (/q that gives 

XQ/27r = 5 MHz, we obtain a Purcell limited qubit life¬ 
time of 2 ms. To next order in 5 q/A the qubit lifetime is 
dependent on the cavity occupation; however, as shown 
in [42], the lifetime increases for higher photon numbers 
in the cavity. 

When the JPM is brought on resonance with the cavity, 
the cavity-JPM states hybridize, so that the eigenstates 
are now 


a) = ;^ i\n, 0) -k (-l)“|n - 1,1)), (C5) 

where the ground state | 0 ,a) = | 0 , 0 ), and the index a 
labels what is normally labelled by ±. In light of this, 
to study qubit relaxation via the JPM, we must examine 
transitions from the two states |l,n,a) to the two states 
|0,n, &) via the JPM-environment coupling operator a^, 
where now 


|l,n,a) 


(C 6 ) 


/ g^in + l)\ 

2A^ jll.n.o)- 


gqVnTT 


A 


jo, n + 1, a), 


\0,n,b) 


(C7) 


1 - | 0 ,n, 6 ) -k ^^^|l,n- 1 , 6 ). 


Similar to the case for cavity-mediated decay, the decay 
rates for these processes are given by 

r^.olA ^ (C8) 

where 7 R(a;) is the JPM’s coupling constant with the 
environment. Using the fact that 

(n-k l,o|cr^|n,6) = n>0 (C9) 

(l,a|aJ|0,6) = ^^, n = 0 (CIO) 

and all other matrix elements are zero, we find that 

(l,n,a|o-^|0,n,6) = {Vn - Vn + l) (Cll) 

= (- 1 )“;^- (C 12 ) 

to first order in ^q/A. As before, assuming an 

Ohmic spectral density we approximate 7 r(^ 077175 ) 

by 7 R(<j-’j)a;Q/a;j, which assumes that the qubit energy- 
shifts due to both the cavity and the JPM are suffi¬ 
ciently smaller than wq (i.e. A^ 3 ^ « ujq). Using 
7 r(|^j) = 200 MHz as in the main text and other quan¬ 
tities as before, we obtain a JPM-limited qubit lifetime 
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« 2 /zs for vacuum in the cavity. For cavity Fock 
states with n > 0 this lifetime scales as 

oc - o = (Vn + Vn + l) , (C13) 

[y/n - y/n + l) 

and so we have oc 2n to leading order. Thus, for 
n = 10 we have an improvement of the qubit lifetime to 
« 40 IIS. 

In addition, there are several competing incoherent 
processes in the JPM, namely bright and dark counts, 
which block the JPM-mediated qubit decay channel. Nu- 


merical simul ations indicate that due to the fact that 
7 j piraa.Onfc relevant photon numbers in the cavity, 
there is almost no appreciable decay of the qubit during 
the measurement process (see Fig. 7(a)). This can be un¬ 
derstood by the fact that oc (/Q/(nA^), which 

for n = 10 is only 0.05%, and so we expect no more than 
a 0.05% change in the qubit state due to JPM-mediated 
decay during the measurement protocol. Detailed study 
of qubit decay during JPM measurement will be the sub¬ 
ject of future work. 
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